Helimagnetism in a disordered and geometrically frustrated quantum spin chain. 
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Neutron diffraction and bulk measurements are used to determine the nature of the low- 
temperature ordered state in a S — 1/2 spin-chain compound with competing interactions. The 
magnetic structure is found to be helimagnetic, with a propagation vector (0.5, £,0), C = 0.174. 
The nearest-neighbor exchange constant and frustration ratio are estimated to be J = 5.8 meV and 
J2/J1 = 0.29, respectively. For idealized spin chains these parameter values would signify a gapped 
spin-liquid ground state. Long-range ordering is attributed to intrinsic non-stoichiometry. 



In recent years the focus in quantum and low- 
dimensional magnetism has shifted from the simplest 
nearest-neighbor chain models to more complex spin net- 
works, effects of impurities, disorder, and inter-chain cou- 
pling. Of particular interest are cases that involve strong 
geometric frustration of interactions. Any of these ef- 
fects may either enhance or suppress quantum fluctua- 
tions, and thus play a crucial role in determining the 
ground state and magnetic properties. Even the simplest 
frustrated models, such as the 5=1/2 spin chain with 
competing nearest-neighbor (nn) next-nearest-neighbor 
(nnn) antiferromagnetic (AF) interactions, have a rich 
phase diagram that includes both gapped and quantum- 
critical gapless phases jjj. Among the latter is, for exam- 
ple, a unique disordered chiral state |2|, which is all that 
quantum fluctuations preserve of helimagnetism, realized 
in the classical version of the model. The challenge is 
to find experimental realizations of frustrated-chain con- 
structs, and to understand the effect of disorder, impuri- 
ties and extended interactions in such materials. 

The charge-ordered compound LiCu2O2 0, 0, S IE 0j 
may, in fact, be just the right model system, featuring 
5=1/2 chains with strong nnn interactions. Paradox- 
ically, LiCu202 possesses characteristics of a quantum- 
disordered "dimer liquid" state, including a gap A sa 
6 meV in the magnetic excitation spectrum |f|, and yet 
undergoes at least one magnetic ordering transition, at 
T c 23 K [1 H 0. The phenomenon is poorly un- 
derstood, and the key missing pieces of information are 
the actual structure of the ordered phase and the mag- 
nitude of frustration. These issues are topics of an ongo- 
ing controversy. /x-SR studies indicate a single complex 
non-collinear state below T c rs 22.5 K, with a precursor 
transition at T± = 24 K. Some recent bulk measurements 
point to an additional transition at T* = 9 K. None 
of the published data provide direct information on the 
frustration ratio of nnn and nn interactions a = Jil J\- 
In the present work we show that in samples with thor- 
oughly characterized stoichiometry the transition at T c 
results in an incommensurate helimagnetic state that per- 
sists to low temperatures. We unambiguously determine 
the frustration ratio, and propose an explanation for the 




FIG. 1: Crystallographic unit cell of LiCu202 showing the 
magnetic Cu 2+ sites (green balls) and the planar helimagnetic 
spin structure (arrows) determined in this work. 



precursor transition at T\. Several alternative mecha- 
nisms of magnetic ordering in what "should" have been 
a quantum-disordered gapped system are proposed. 

As discussed in detail in Refs. 0, 0, 0] , LiCu202 con- 
tains an equal number of Cu + and Cu 2+ ions in distinct 
non-equivalent crystallographic positions. The magnetic 
Cu 2+ ions carry 5 = 1/2 and form "triangular" two- 
leg ladders (Fig. [TJ, which can also be viewed as zig-zag 
chains with competing nn and nnn interactions, J nn and 
J nnn , respectively. These chains run along the b axis of 
the orthorhombic crystal structure, and are well sepa- 
rated from each other by double chains of non-magnetic 
Li + ions and layers of non-magnetic Cu + sites. One key 
element of the present work was the preparation of sam- 
ples with a thorougly controlled chemical composition. 
Single crystals of LiCu202were grown in an alundum cru- 
cible in air atmosphere using the self- flux method. The 
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FIG. 2: Temperature dependence of magnetic susceptibil- 
ity of LiCu2C>2 measured in a magnetic field H — 100 Oe 
applied parallel to the (a, b) cleavage plane (symbols). The 
solid line is a fit based on the frustrated S — 1/2 chain model, 
as described in the text. Taking a numerical derivative (inset) 
reveals a phase transition T c — 22 K (arrow). 

lattice parameters a = 5.730(1)A, b = 2.8606(4)A, and 
c = 12.417(2)1 were verified by powder X-rays diffrac- 
tion, which also confirmed the absence of any apprecia- 
ble amounts of impurity phases. These values are in very 
good agreement with those reported in All crystals 
were found to be microscopically twinned with respect to 
the [1,1,0] plane, so that a rj 26. 

Unexpectedly, it was found that those samples with 
the smallest mosaic spread and most reproducible mag- 
netic behavior had lower content of Cu + ions than fol- 
lows from the stoichiometic formula Lii+Cu2+C>2. The 
density was determined to be p — 5.12 g/cm 3 at room 
temperature, which corresponds to an actual composition 
Li1.i6Cu1.g4O2.01- Chemical disorder and a Cu-deficiency 
by as much as x — 16 % are thus inherently present. 
The "surplus" Li + ions in Li1.i6Cu1.s4O2. 01 occupy Cu 2+ 
sites, due to a good match of ionic radii. A substitu- 
tion of Cu + ions by Li + ions in the LiCu2 02 structure 
is prevented by the dumbbell oxygen coordination be- 
ing very uncharacteristic for the latter ion type. Charge 
compensation requires that the introduction of 16% non- 
magnetic Li + ions into the double chains is accompa- 
nied by a transfer of 16% of the S = 1/2-carrying Cu 2+ 
ions onto the Cu + inter-chain sites. As will become cru- 
cial for the discussion below, our disordered and non- 
stoichiometric Li1.i6Cu1.s4O2. 01 crystals (referred to as 
simply "LiCu202" throughout the rest of the paper) have 
appreciable concentrations of both non-magnetic Li + im- 
purities in the zig-zag chains, and magnetic Cu 2+ impu- 
rities positioned in-between chains. 

The single crystal samples were characterized by bulk 
susceptibility and specific heat measurements. x(T) 
data were taken in a commercial SQUID magnetome- 
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FIG. 3: The measured temperature dependence of specific 
heat in LiCu2 02(open symbols) indicates a phase transition 
at T c ~ 22 K and a possible precursor at Ti w 25 K (arrows). 
Subtracting the phonon contribution (solid line) allows to ex- 
tract the temperature dependence of magnetic entropy (lower 
right inset). 



tcr in the temperature range 5-100 K and a magnetic 
field H = 100 Oe applied parallel to the (a, 6) cleav- 
age plane and are shown in symbols in Fig. [21 The 
main feature is a broad maximum at T w 36 K char- 
acteristic of a quasi-one-dimensional magnet, that signi- 
fies the formation of short-range correlations within the 
chains. Taking the temperature derivative of the mag- 
netic susceptibility (Fig. [21 inset) reveals a sharp anomaly 
at T c = 22 K, that we attribute to the onset of long- 
range magnetic order. The high-temperature part of the 
x(T) curve is expected to be uninfluenced by this phase 
transition, representing the behavior of individual zig- 
zag chains. The data taken above T = 50 K were ana- 
lyzed using the high-temperature expansion formulas for 
the susceptibility of the quantum S = 1/2 frustrated- 
chain model 0. An excellent fit (solid line in Fig. |2J) 
is obtained with J\ = 5.9(1) mcV, g = 2.114(2), and 
a = 0.282(3). The anomaly at T c is also manifest in the 
specific heat data measured using a "Termis" quasiadia- 
batic microcalorimeter and plotted in Fig. [31 The peak 
observed at T « T c is well-defined, but, as indicated 
by arrows in the blowup plot, actually has a charac- 
teristic flat top that extends between T c = 22 K and 
Ti = 24 K, in agreement with the results of Ref. 0. 
The solid line in Fig. [31 represents a crude estimate for 
the phonon contribution, calculated under the assump- 
tion that the total magnetic entropy released at short- 
range and long-range ordering in LiCu202 is equal to 
ii In (25 + 1) = 5.76 J/molc K. The temperature depen- 
dence of magnetic entropy extracted in this fashion in 
plotted in the lower right inset of Fig.[3J A large fraction 
of the entropy is released above the ordering temperature, 



as expected for a low-dimensional system. The sharp 9 K 
anomaly reported in Ref. is totally absent in our x(T) 
and C(T) data. We suspect that this feature is due to an 
impurity phase, most likely Li 2 Cu02, which is known to 
go through an AF transition at 9 K . Otherwise, the 
characteristics of our samples seem to be in good agree- 
ment with those reported by other authors 0, S IS • 

The nature of the magnetically ordered state was de- 
termined in a single-crystal neutron diffraction experi- 
ment. Since LiCu202 was prepared using a naturally oc- 
curing Li isotope mixture, the attenuation of a neutron 
beam due to absorbtion by 6 Li nuclei was significant. 
The sample was cut parallel to the (a, b) crystallographic 
plane to the shape of a thin plate 0.9 x 15 x 15 mm 3 . The 
measurements were performed at the HB1 and HB1A 
3-axis spectrometers installed at the High Flux Isotope 
Reactor at ORNL. For measuring integrated Bragg in- 
tensities the instruments were used in 2-axis mode, with 
a well-collimated incident neutron beam of a fixed en- 
ergy Ei = 14.7 mcV. A Pyrolitic Graphite PG(002) re- 
flection was used in the monochromator. A 5 cm thick 
PG filter was positioned in front of the sample to elim- 
inate higher-order beam contamination. Bragg intensi- 
ties were collected in rocking curves and corrected for 
the usual Lorentz factor. Absorbtion corrections were 
applied assuming a thin-plate geometry. High-resolution 
measurements of the magnetic propagation vector were 
performed in 3-axis mode, with a PG(002) analyzer and 
48' - 40' - 40' - 240' collimators. 

The main finding of this work is that below T c w 22 K 
LiCu?Q? acquires incommensurate magnetic long range 
order 18] . The phase transition leads to the appear- 
ance of new Bragg reflections that can be indexed as 
(2^1, fe ± 77, 1), n, k, Wnteger, ( « 0.174. Such peaks 
were observed in both crystallographic twins. The peak 
widths were found to be resolution- limited along all three 
crystallographic directions at all temperatures below T c . 
A typical scan across the (0.5, 0.826, 0) reflection taken at 
T = 2 K is shown in the inset of Fig.0k. The correspond- 
ing peak intensity is plotted as a function of temperature 
in Fig.^Ji. A simple power law fit to the data taken above 
T = 15 K yields T c = 22.3(6) K and (3 — 0.25(0.07). The 
residual intensity seen in Fig. at T > T c is due to 
critical scattering and is quite broad in q-space. 

The temperature dependence of the magnetic propa- 
gation vector was deduced from Gaussian fits to /c-scans 
across the (0.5,0.826,0) peak, and is plotted in symbols 
in Fig. 0Jd. Below T ks 17 K the incommensurability 
parameter is practically T-independent and appears to 
have a strictly incommensurate value ( = 0.1738(2). In- 
terestingly, as T c is approached from below, Q progres- 
sively decreases as indicated by the arrow in Fig. 
The minimum value of Q observed in our experiments is 
about 0.172. The temperature dependence of C, may be 
the cause for a flat-top peak in the measured specific heat 
curve. In the close vicinity of T c , where £ changes rapidly, 



>, 
'(0 



o 



ID 
O 



LiCu O 

2 2 



T =22.3(0.6) K 




0.174 



0.172 



(b) 




/X 





10 15 20 

T(K) 



25 30 



FIG. 4: (a) Measured temperature dependence of the 
(0.5, 0.826, 0) magnetic peak intensity in LiCu202 (symbols). 
The solid line is a power law fit to the data. Inset: fc-scan 
across this reflection measured at T — 2 K. (b) Measured in- 
commensurability parameter £ plotted as a function of tem- 
perature. The solid line is a guide for the eye. 



the large-period structure may undergo one or more 
"devil's staircase" -type commensuration transitions, par- 
ticularly at Ti 0- Intensity being the limiting factor, we 
were unable to locate any well-defined magnetic Bragg 
relections in the narrow temperature range T c < T < T\ . 

The spin arrangement in the ordered state was de- 
duced from the analysis of 23 non-equivalent magnetic 
Bragg peaks with < h < 3.5, < k < 1.5 and 
< / < 8. Those reflections for which the incident or 
scattered beams formed angles of less that 15° with the 
crystalline plate were not included in this data set, since 
the corresponding absorbtion corrections are expected to 
deviate substantially from the thin-plate approximation. 
After examining a number of models, we found that ex- 
cellent fits to the data are obtained assuming a planar 
spin- helix propagating along the zig-zag Cu 2+ chains, 
with a fixed relative rotation angle <j) = tt(1 — £) between 
consecutive spins. Any spins related by a translation 
along the c axis and b axes are parallel and antiparallcl 
to each other, respectively. The magnetic structure factor 
for such a helical state can be easily calculated analyti- 
cally [ill ]. After averaging over possible g-domains and 
taking into account the magnetic form factor of Cu 2+ , 
an excellent fit to the measured magnetic intensities is 
obtained assuming all spins to be confined to the (a, b) 
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crystallographic plane, and adjusting only an overall in- 
tensity scaling factor. This spin structure is visualized 
in Fig. ^ Unfortunately an independent estimate of the 
magnitude of the ordered moment from our data was pre- 
vented by rather severe extinction corrections of nuclear 
Bragg intensities typically used for normalization. 

A helimagnetic state is the typical way geometric frus- 
tration is resolved in a classical magnet. Perhaps the 
first of many known examples is that of M11O2, dat- 
ing back to the 50's 0]. The incommesurability pa- 
rameter £ is directly related to the frustration ratio a 
through Act = 1/ cos(7r£). In our case, for £ = 0.174 
this yields a = 0.29, in excellent agreement with our 
estimates based on high-temperature bulk susceptibil- 
ity data. What makes LiCu202 remarkable, is that ac- 
cording to existing theories for quantum S — 1/2 spin 
chains with nnn interactions, this value corresponds to 
a gapped disordered dimer-liquid state, not too far from 
the collinear gapless phase realized for a < 0.241 |l3|. 
This would seem to agree well with the observation of a 
spin gap at T > T c 6] , but is of course incompatible with 
incommensurate long-range order at lower temperatures. 
One possible explanation for the observed behavior are 
strong inter-chain interactions, such as those that induce 
ordering in the well understood Haldane-gap compound 
CsNiCl3 0] . Here the amplitude of transverse dispersion 
of the energy gap increases with decreasing T, eventually 
driving the gap to zero at the 3D zone-center and induc- 
ing a soft-mode phase transition. A similar mechanism 
may be active in LiCu2 02, and is not incompatible with 
the ESR data of Ref. |6j. Nevertheless, it is difficult to 
envision the necessary strong inter-chain superexchangc 
pathways in the stoichiometric material. 

Chemical disorder, that is inherently present in our 
samples may be expected to play an important role in 
magnetic ordering. In particular, the 16% Li + non- 
magnetic defects in the double chains will locally de- 
stroy valence bond spin-singlets and liberate end-chain 
spin degrees of freedom. The free spins can then order 
in three dimensions via arbitrary small inter-chain in- 
teractions. This mechanism is known to drive mag netic 
ordering in the spin-Peierls compound CuGe03 |lj| , and 
the Haldane-gap system PbNi 2 V20s 0. Alternatively, 
long-range ordering in LiCu2 02 may be driven by the 
16% Cu 2+ impurities positioned in-between the gapped 
spin chains. Such magnetic impurities carry interactions 
across the non-magnetic Cu + layers, completing a 3D 
network of magnetic interactions. Long range order oc- 
curs for an arbitrary small concentration of impurities, 
and involves ordering of both the impurity spins and the 
chain spins. This model in realized in (^Yi-x^BaNiOs 
rare earth nickelates, where classical R 3+ ions bridge the 
Ni-based gapped spin chains 0]. To fully understand 
the role of impurities in the ordering of LiCu202 more 
detailed neutron scattering experiments on isotopically 
pure samples are required. 



In summary, we have observed incommensurate heli- 
magnetism in a S — 1/2 zig-zag chain compound where 
the estimated frustration ratio corresponds to gapped 
disordered phase. We speculate that long-range magnetic 
ordering is aided by intrinsic impurities and deviations 
from stoichiomctry. This work was partially supported 
by RFBR Grants 02-02-17798 and 03-02-16108. Work at 
ORNL was carried out under DOE Contract No. DE- 
AC05-00OR22725. 
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